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ABSTRACT 

We use axisymmetric magnetohydrodynamic (MHD) simulations to investigate the 
launching and collimation of jets emerging from the disk-magnetosphere boundary 
of accreting magnetized stars. Our analysis shows that the matter flows into the jet 
from the inner edge of the accretion disk. It is magnetically accelerated along field 
lines extending up from the disk and simultaneously collimated by the magnetic pinch 
force. In the reference run which we use for analysis, the matter in the jet crosses 
the Alfven surface a few i?* above the disk and the fast magnetosonic surface ~13 
above the disk. At larger distances, the magnetic pressure is a few times smaller 
than the total matter pressure but the magnetic force continues to accelerate and 
collimate the jet. In steady state, we observe a matter ejection-to-accretion ratio of 
~0.2. Across different simulation runs, we measure a range of half-opening angles 
between 6 w 4° and 20° at the top of the simulation region, depending on the degree 
of magnetization in the outflow. We consider the case of stars undergoing epochs of 
high accretion (such as EXors, FUORs, and CTTSs) where the stellar magnetosphere 
is strongly compressed by the incoming accretion disk. For a typical EXor (mass 0.8 
Mq, radius 2 R©) accreting at 10~^ MQ/yr, we measure poloidal velocities in the jet 
ranging from 30 km/s on the outer edge of the jet to more than 260 km/s on the inner 
edge. In general, the models can be applied to a variety of magnetized stars — white 
dwarfs, neutron stars, and brown dwarfs — which exhibit periods of high accretion. 
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1 INTRODUCTION 

Young, accreting stars produce strongly collimated high-speed 
outflows (jets) which are a key mechanism in transporting 
mass, energy and angular momentum out of the disk and 
facilitating accretion onto the young protostellar core. The 
formation of a collimated jet requires tandem processes: the 
matter must be accelerated out of the disk and undergo col- 
limation into a jet. A number of theoretical models have 
been proposed to explain the mechanism which launches the 
matter from the disk (see review by Ferreira et al. 2006). 
The first type is the magneto-centrifugal mechanism in which 
the outflow is launched from the disk as an extended wind 
(Blandford & Payne 1982; Konigl & Pudritz 2000). A second 
class of models involves outflows driven from the the disk- 
magnetosphere boundary, either as an X-wind type outflow 
(Najita & Shu 1994; Shu et al. 1994; Cai et al. 2008) or as a 
magnetically launched conical wind (Romanova et al. 2009, 
hereafter referred to as R09). The matter may also be magnet- 
ically launched from the inner disk (Lovelace et al. 1991) or 
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driven as a stellar wind (Matt & Pudritz 2008). Observations 
of stars with strong outflows show that the typical outflow ve- 
locities are of the order of the Keplerian velocity of the inner 
disk region, favoring the models where the outflow originates 
from the disk-magnetosphere boundary or the inner disk. As 
the matter flows out of the disk, helical magnetic field lines 
frozen into the rotating outflow can self-coUimate the mat- 
ter into a jet (Lovelace, Wang & Sulkanen 1987). The out- 
flow may also become collimated by interaction with an exter- 
nal poloidal field threading the disk (Matt et al. 2003; Fendt 
2009). Yet another possibility is that the matter becomes col- 
limated by pressure from an external medium (Lovelace, Berk 
& Contopoulos 1991; Frank & Mellema 1996). 

Spectral measurements show that a significant number 
of Classical T Tauri Stars (CTTSs) exhibit signs of outflows 
from their disks (Edwards et al. 2003, 2006; Gomez de Cas- 
tro & von Rekowski 2010) with early observations indicating 
a correlation between the outflow and mass accretion rates 
of CTTSs (Cabrit et al. 1990; Hartigan et al. 1995). Even 
with improvements to telescopes, direct imaging of the in- 
ner disk region is still hampered by the insufficient resolution 
of ground- and space-based observatories (see review by Ray 
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et al. 2007) and observations have only recently shown that 
the jets become coUimated at distances less than about 10 AU 
from the star (Hartigan et al. 2004; Coffey et al. 2008). In lieu 
of direct imaging of the inner disk, astronomers have relied on 
numerical simulations to study the structure and dynamics of 
the disk-magnetosphere boundary. 

Simulations Early simulations performed by Hayashi, 

Shibata & Matsumoto (1996) and Miller & Stone 
(1997) achieved single-episode outflows from the disk- 
magnetosphere boundary lasting a few dynamical timescales. 
As simulations grew more sophisticated, longer runs per- 
formed by Goodson et al. (1997, 1999), Hirose et al. (1997), 
Matt et al. (2002), and Kiiker et al. (2003) showed several 
episodes of field line inflation and outflows. However, none 
of these early simulations produced robust outflows which 
lasted long enough to estabUsh the outflow behavior and 
dynamics. Long lasting outflows have been achieved by 
treating the disk as a boundary condition and launching 
matter into the corona (e.g. Romanova et al. 1997; Ouyed 
& Pudritz 1997; Ustyugova et al. 1999; Krasnopolsky et al. 
1999; Fendt & Elstner 2000; Matsakos et al. 2008; Fendt 
2009; Staff et al. 2010). However, In order to understand the 
launching and coUimation mechanisms of the outflows, the 
simulations must include a realistic, low-temperature accre- 
tion disk and solve the full magnetohydrodynamic (MHD) 
equations in both the disk and coronal space. Recently, there 
has been much work in this direction with the modeling of 
outflows launched from realistic disks threaded by large scale 
magnetic fields (Casse & Keppens 2002, 2004; Zanni et al. 
2007; Murphy et al. 2010). 

Our group first obtained long lasting outflows from re- 
alistic accretion disks for the case of rapidly rotating stars 
in the "propeller" regime (Romanova et al. 2005; Ustyugova 
et al. 2006) . These simulations showed a two-component out- 
flow in which most of the matter is carried away through a 
conical-shaped wind while most of the energy and angular 
momentum flows into a low-density, high-velocity axial jet. 
This axial jet is magnetically-dominated and well coUimated 
by the toroidal magnetic field. More recently, we observed 
single-component, long lasting conical outflows in the more 
general case of slowly rotating stars (R09). In R09, the coni- 
cal outflows are driven from the disk-magnetosphere bound- 
ary by magnetic pressure and unlike X-winds, these outflows 
do not require equality of the magnetospheric and corotation 
radii and can originate even from slowly rotating stars. How- 
ever, the conical winds showed only weak coUimation within 
the simulation region. 

In this work, we present axisymmetric MHD simu- 
lations of long-lasting, coUimated jets launched from the 
disk-magnetosphere boundary of a magnetized star with an 
aligned dipole. Building upon the previous work on conical 
winds by R09, we simulate an accretion disk with a higher 
accretion rate in a larger simulation region and observe the 
emergence of a magnetically driven, magnetically coUimated 
jet. In §2 we describe the numerical methods and initial and 
boundary conditions used in the simulations. Then, in §3, we 
give an overview of the simulations and describe the reference 
run which we use for analysis of the launching and coUima- 
tion mechanisms. In §4 we describe the fluxes, velocities and 
forces in the jet and investigate the jet launching and coUi- 
mation mechanisms. Lastly, in §5 we discuss the relation of 



this work to previous results as well as the application of our 
simulations to several tjqses of rapidly accreting young stars. 



2 NUMERICAL SETUP 

Here, we briefly describe the main aspects of our numeri- 
cal model. The numerical model used here is identical to the 
setup used in R09 and utilizes an axisymmetric Godunov-type 
code based on the equations of viscous, resistive MHD. In the 
inertial reference frame, these equations are 



dt 

d{pv) 

dt 



V • (pv) = 0, 

+ v-r = pg, 



dB 

dt 



-Vx(vxB) + Vx (r7tV x B) = 0, 



djpS) 
dt 



+ V • {pSv) = Q. 



(1) 
(2) 
(3) 
(4) 



In this set of equations, p is the density, v is the flow velocity, 
T is the momentum flux-density tensor, g = —{GM/r^)r is 
the gravitational acceleration due to the star, B is the mag- 
netic field, rjt is the magnetic diffusivity coefficient (see Ap- 
pendix B), S is the specific entropy, and Q is the rate of 
change of entropy per unit volume. We treat the plasma as 
a monatomic ideal gas such that S = Cv ln(P/p'') where P 
is the gas pressure and the adiabatic index 7 = 5/3. 

Grid The code uses a spherical coordinate system (r, 6*, 0) 
where is the colatitude and 4) is the azimuth angle. Ax- 
isymmetry imposes the additional condition d/d(b = on the 
equations of MHD^. Since we use a Godunov scheme, all nu- 
merically calculated variables are cell centered on the grid, 
except for the magnetic vector potential which is calculated 
on the nodes. The simulations are performed in the region 
Rt < r < Rout, < 9 < Tv/2 and reflected across the 
axisymmetry axis for plotting (the "z-axis" in cylindrical co- 
ordinates). The grid is uniform in the 9 direction and the 
size steps in the radial direction are chosen such that the 
poloidal-plane cells are curvilinear rectangles with approxi- 
mately equal lengths on each side. This choice results in high 
spatial resolution near the star where the disk-magnetosphere 
interaction takes place while also permitting a large simula- 
tion region. In the simulations presented in this paper, we 
choose the number of grids in the 9 and r directions to be 
No = 50 and A''^ = 120, respectively. This corresponds to 
d6 = n/Ne = 1.8° and a total region radius of about 42 R*. 
The smallest grid size at the surface of the star is (7r/100)R* 
on a side. 



2.1 Dimensionless variables 

Within the numerical code, the equations of MHD are 
reparametrized with normalized variables (e.g. p = p/po, 
B = B/Bo, V = v/vq, etc.) and solved in a dimension- 
less form. To further simplify the equations, we also take 
GM« = 1 and 71=1. The dimensionless equations permit 



^ The equations of axisymmetric MHD in spherical coordinates can 
be found in Ustyugova et al. (2006). 
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us to apply the general results to a wide variety of accreting 
stars. The calculation of the reference units is given in Ap- 
pendix A and sample reference values for EX Lupi (EXor), FU 
Orionis (FUOR) and CTTS class stars are shown in Tab. Al. In 
order to convert the dimensionless units into real values, mul- 
tiply the dimensionless value by the corresponding reference 
value in Tab. Al. For example, to calculate t = 1000 for a typ- 
ical CTTS, multiply by the reference time (to.CTxs = 0.366 
days) to get t = 366 days. 

For the remainder of the paper, all values and variables 
are given in terms of the normalized units (with tildes im- 
plicit) except where explicitly assigned physical units. 



2.2 Numerical method 

To numerically integrate the MHD equations, we split the 
physical processes into four blocks: (1) an "ideal MHD" 
block in which we calculate the dynamics of the plasma and 
magnetic field without dissipative processes; two "diffusion" 
blocks (2) and (3) in which we calculate the diffusion of the 
poloidal and azimuthal components of the magnetic field for 
frozen values of the plasma velocity and thermodynamic pa- 
rameters (density and pressure); and a "viscosity" block (4) 
in which we calculate the viscous dissipation due to the rep 
and 64) components of the viscous stress tensor^ . Integration 
of the equations in time is performed with a two-step Runge- 
Kutta method. To determine the fluxes between the cells, we 
use an approximate solution of the Riemann problem anal- 
ogous to the Roe solver described in Brio & Wu (1988), ex- 
cept that we take the energy conservation equation in the en- 
tropy form. The dynamical variables are determined in the 
cells while the vector-potential of the magnetic field, A^, is 
determined on the corner nodes. We guarantee the absence 
of magnetic charge by calculating Aj, at each time step and 
then using it to obtain the poloidal components of the mag- 
netic field (Br, Be) in a divergence-free form (Toth 2000). 
This ensures that the divergence-free condition V • S = is 
always satisfied to within machine accuracy everywhere in the 
simulation region. We set a floor density of pfir — 2.5 x 10~^ 
throughout the simulation region to prevent the density from 
vanishing near the axis. While the jet is being launched, this 
adds a small amount of matter to the grids nearest the axis. 

For a complete description of the numerical method as 
well as standard tests of the code, see Appendices C and D in 
R09. 



2.3 Initial and boundary conditions 

Initial Conditions We take the stellar magnetic field to be 
an aligned dipole B, = [3(/Lt • r)r — p.r^]/r^ where is 
the stellar dipole moment which is taken as a parameter in 

the simulations. Initially, the entire simulation region is filled 
with a non-rotating low-density, high-temperature isothermal 
plasma and there is no disk present. In spherical coordinates, 
the initial density and pressure distributions are 

p = Pc exp[GM/(7^Ter)], P = exp[GM/(7^^er)], (5) 



^ See Appendix B for an in-depth discussion of the treatment of the 
turbulent viscosity and diffusivlty in the simulations. 



where Tc is the coronal temperature and pc is the coronal den- 
sity at the external boundary, related by the ideal gas equa- 
tion of state Pc = pcTZTc. The star initially corotates with the 
outer boundary such that = {GM^/Rt^^Y^'^. Between t 
= and t = 100 the star is gradually spun up, reaching a final 
value of 0,j = [GM./(37?.)^]^/^ = 0.19 corresponding to a 
disk corotation radius Rcot = 3. Information about the stellar 
rotation propagates out at the Alfven speed along the stellar 
magnetic field lines into the low-density corona. 

External boundary The external boundary i?out is divided 
into a coronal region (0 < < 9d) and a disk region id^ < O 
< it/2) where, again, is the colatitude measured from the 
axisymmetry axis. The disk height at the boundary is set by 
satisfying hydrostatic equilibrium. Starting at t = 0, we permit 
low-temperature, high-density matter to flow in through the 
disk boundary region (0 > 6d) with a fixed density profile set 
by pd- The matter velocity along the disk boundary is set to 
be slightly sub-Keplerian so that the matter will flow into the 
simulation region. All other hydrodynamic variables have the 
"free" boundary conditions in the disk region, d{. . .)/dr = 0. 

In the coronal region (6 < da), the boundary conditions 
are also free for all hydrodynamic variables. However, we pro- 
hibit matter from flowing into the simulation region from this 
portion of the boundary. We solve the transport equation for 
the flux function $ so that the magnetic flux flows out of the 
region together with matter. 

The boundary conditions on the equatorial plane and 
on the rotation axis are symmetric and axisymmetric, respec- 
tively. On the axisymmetry axis, we enforce the axisjmimetry 
boundary conditions by requiring that the toroidal variables 
obey B4,{-9) = -84,(6), V4,{-9) = -V4,{e), et cetera. This 
imposes the explicit requirement that the toroidal variables 
must cross on the axis. Nonetheless, we must note that there 
is a region near the axisymmetry axis (9 = to 6c) where the 
azimuthal angular velocity v^i/r may take on spurious values. 
This can in turn result in spurious values of near the axis 
even though 8^,(6 = 0) = 0. However, as we will discuss in 
§3.2, the anomalous values do not affect the collimation of 
the jet. In our simulations 6c spans about two angular grid 
cells or about 3.6°. 

Internal boundary on the star The inner boundary Ri^ = 
R* lies on the stellar surface. We assume that the poloidal 

component of the magnetic field Bp is frozen to the surface 
of the star: Br is held fixed while Bg and obey the free 
boundary conditions dBe/dr = and dB^^/dr = 0. The 
plasma pressure and entropy also obey free boundary con- 
ditions on the stellar surface; at each timestep, the density on 
the surface is recalculated from these two free variables. The 
matter velocity components are calculated using free bound- 
ary conditions and then adjusted to be parallel to the mag- 
netic field vectors on the stellar surface. We do not consider 
outflows due to stellar winds and therefore we only permit 
matter to flow inward onto the star. 



3 SIMULATION OVERVIEW 

In order to investigate formation of jets from the disk- 
magnetosphere boundary, we performed dozens of simula- 
tions with varying initial parameters and grid sizes. The re- 
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Figure 1. The jet at time t = 860. The background shows the poloidal matter flux density p | Vp | and the lines show contours of the magnetic flux 
function which is a proxy for the poloidal magnetic field. The horizontal axis shows the distance away from the star in units of the reference 
radii Rq . The red vectors show the poloidal matter velocity Vp . To obtain dimensional values, multiply these numbers by the numbers given in 
Tab. Al. For a CTTS, to = 0.356 days and Rg = 2RQ(e.g. column 1 in Tab. Al): in real units then, the time t=860 corresponds to 315 days and 
the simulation region is 0.39 AU in radius. 
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Figure 2. Formation and coUimation of the jet at different times t. The p |vp| background contour and field lines are plotted on the same scale as 
the contours in Fig. 1 . 



gion size is set by Ng and Nr, the number of grid cells in the 
9 and r directions, respectively. The simulations have a total 
of 5 free parameters: and ad, the dimensionless viscosity 
and diffusivity coefficients; pd and pc, the boundary densi- 
ties of the disk and corona; and fi, the stellar magnetic dipole 
moment. A sample of the runs which successfully produced 
strongly coUimated jet-like outflows is shown in Tab. 1. 

We use the Keplerian rotation period at the surface of 
the star to as a reference timescale. We observe three general 



classes of outflows: one-time episodic outflows, stable weakly 
collimated conical-type outflows, and stable collimated jets. 
The one-time episodic outflows occur when the accreting mat- 
ter strongly compresses the stellar field lines, driving transient 
outflows from the disk-magnetosphere boundary. In this case, 
the matter and magnetic pressures quickly re-equilibrate, 
halting the driving mechanism and suppressing the outflow 
in just a few to- The conical winds are similar to those previ- 
ously studied in R09 and are driven by the compressed mag- 
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Table 1. Parameters for simulation runs which produced well coUi- 
mated jets. Ng and Nr are the number of grids in the 9 and radial di- 
rections, respectively; a„ and are the dimensionless viscosity and 
diffusivity coefficients. In dimensionless units, is the fixed density 
on the disk boundary; pc is the fixed coronal density at the external 
boundary; -Rout is the radius of the outer simulation boundary (this 
is set by Ng and A''^, so it is not a free parameter), and p is the stellar 
magnetic dipole moment. The bold line is the reference case which 
we use for analysis in this paper 

netic field at the disk-magnetosphere boundary. In this work, 
we are interested in particular in the last class of observed 
outflow: the stable collimated jet. 

3.1 Reference run 

Fig. 1 shows one of the simulations which produced a strongly 
collimated jet: throughout the remainder of the paper, we will 
use this run as a reference case to perform additional analy- 
sis. The background contours represent the poloidal matter 
flux p|vp| and the lines represent contours of the magnetic 
flux function >]/, which obeys the relation B • Vip — 0. In 
the poloidal plane, therefore, the ip contours represent the 
magnetic field lines. The reference simulation uses a spheri- 
cal grid with Ng = 50 cells in the 6 direction and A'^,. = 120 
cells in the radial direction with the parameters: = 0.3, 
Od = 0.1, disk density pd = 5, and coronal density pc — 0.001 
(a summary of the parameters is shown as the bolded line in 
Tab. 1). For analysis of the launching and collimation mecha- 
nisms, we study the simulation at t = 860 when the jet is well 
established. 

Fig. 2 shows snapshots of the reference simulation at six 
moments in time ^. At t = 100, the matter flowing in from the 
boundary begins compressing and inflating the dipole field 
around the star. By t = 300, the disk is only a few stellar 

^ An animation of the full reference simulation is available 
at http://astro.cornell.edu/~pslii/research.htm. We encourage the 
reader to refer to the online animations for clarification of the dis- 
cussion in this section. 




Figure 4. 3D view of the poloidal matter flux density p \vp\ contours 
of the jet at t = 860. Three matter flux density surfaces are plot- 
ted: p \vp\ = 0.001, 0.003, and 0.01. The red streamtraces show the 
strongly wound magnetic field lines in the corona of the star which 
coUimate the outflowing matter 

radii away from the star; the close proximity of the disk per- 
mits matter to accrete directly onto the star by flowing along 
the closed dipole field lines. At the same time, the inflated 
field lines in the corona undergo forced reconnection, eject- 
ing plasmoids into the stellar corona. 

The strong magnetic gradient at the disk-magnetosphere 
boundary launches matter along the inflated dipole field lines. 
The outflow first emerges from the open field lines at the 
disk-magnetosphere boundary starting around t = 350 and 
becomes stable and well established by t = 500. As the jet 
stabilizes, the magnetic reconnection within the jet halts and 
the matter flows smoothly. Fig. 3 shows a zoomed-in view 
of the outflow at t = 860, illustrating the proximity of the 
launching region to the star. As the outflow becomes stronger, 
we observe a gradual magnetic collimation of the outflowing 
matter into a jet (see frames t = 700 to t = 1100 of Fig. 2). 

In Fig. 4, we show a three-dimensional rendering of 
the jet at t = 860. Three poloidal matter flux surfaces are 
shown along with tightly wound helical magnetic field lines — 
represented by the red ribbons — which play a key role in col- 
limating the outflow (we discuss this later in §4.3). The jet 
has the shape of an inverted tapered cone with a cylindrical 
radius of 1 — 4 i?, at the base of the jet and 3 — 8 i?, at the 
top. For this reference case, we measure a half-opening angle 
of 6 PS 4° at the top of the simulation region. 

3.2 Properties of the jet 

In order to describe the relative magnitudes of the matter 
and magnetic pressures around the star, we use the standard 
plasma parameter /? and the kinetic plasma parameter /?i : 

The kinetic plasma parameter f3i takes into account the ram 
pressure of the gas, pv^ in addition to the thermal gas pres- 
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Figure 3. A zoomed-in plot of the inner disk region at the time t=860. The plotted contours show the same scales as in Fig. 1. 




5 10 15 20 25 30 35 40 

r 

Figure 5. p |vp| contours overplotted with the /3 = 1 and /3i = 1 lines 
(see Eq. 5). The region interior to the /3i line is matter dominated. 
The dashed white line shows the radius where the 9 cross-sections 
shown in Fig. 5 are taken. Also plotted are a few of the velocity 
streamlines in the jet. 

sure, P. Regions with /?, /^i » 1 are matter pressure dom- 
inated while regions with /3, /3i <C 1 are magnetic pressure 
dominated. Fig. 5 shows the /?, /3i = 1 lines for the t = 860 
reference case. The conventional plasma parameter (3 is much 



less than unity everywhere except within the disk; hence if 
we use the plasma /3 criterion, we find that the jet is com- 
pletely magnetically dominated. However, if the ram pressure 
is taken into account and /3i is used instead, then we find that 
/3i > 1 inside the jet: in other words, the magnetic pressure is 
only a few times smaller than (or comparable to) the effective 
matter pressure, showing that the magnetic field is important 
in driving and coUimating the jet even at large distances from 
the star. 

Fig. 6 shows cross-sections of the jet in the 9 direction 
for angles 6 = 0°-45° at the r = 20 surface (see Fig. 5). The 
dashed vertical lines represent the /3i = 1 surfaces which very 
approximately delimit the boundaries of the jet. We indicate 
the region where numerical artifacts may play a role by grey- 
ing out the two grids closest to the axis. Fig. 6f shows the /3 
and/3i profiles of the jet, with /3 10"" -10"^ and /?i ^ 1-6 
inside, indicating that the jet is "weakly" matter dominated. 
Fig. 6a shows the velocity cross-section of the jet: within the 
jet, the matter is supersonic and the poloidal matter velocity 
Vp is larger than the poloidal Alfvcn speed VAp = | Bp | / y/4¥p, 
except at the left edge of the jet where they are comparable. 
The poloidal velocities in the jet range from Vp ^ 1 on the 
inner edge of the jet to Vp 0.1 on the outer edge. 

Panel 6a also shows that the toroidal matter velocity is 
very large even outside the jet, contributing to the winding of 
the magnetic field lines in the corona. Fig. 6b shows the pro- 
file of the normalized poloidal current Ip = r sin 9B^ flowing 
out of the simulation region through the jet (see Appendix 
C for analysis of the current flow). Figs. 6c and 6d show the 
matter and angular momentum flux profiles in the jet: we see 
that the angular momentum transport is dominated by the 
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Figure 6. Cross-sections of the jet in the 9 direction at r = 20 and t = 860 (see Fig. 5). The dashed vertical hnes indicate the /3i =1 surfaces 
which approximately represent the boundaries of the jet. We indicate the region where numerical artifacts may play a role by greying out the two 
grids closest to the axis. Panel (a) shows the velocity cross-section in the jet. and Vp are the velocities in the toroidal and poloidal directions, 
respectively; Veac is the escape velocity at r = 20; Cs is the sound speed; and VAp is the poloidal Alfven speed. Panel (b) shows the poloidal 
current Ip in the jet. Panel (c) shows the density p and matter flux p | Vp | profiles in the jet. Panel (d) shows and Lf, the angular momentum 
fluxes from matter and magnetic fields, respectively. Panel (e) plots the magnetic field components in the poloidal Bp and toroidal directions. 
And lastly. Panel (f) shows the plasma parameters in the jet (described in Eq. 6). 



outflowing matter* (this is further discussed in §4.1). Panel 6e 
shows the magnetic field distribution; note that the toroidal 

field is much larger than the poloidal field everywhere except 
on the inside edge of the jet, where we observe spurious val- 
ues of B^. Note, however, that the magnetic collimation force 
is proportional to the transverse gradient of Ip = (r sin OB^)"^ 
which is well behaved near the symmetry axis, as shown in 
Fig. 6b, and not on B^ itself (Lovelace et al. 1989). There- 
fore, the spurious axis values of B^i, do not affect the collima- 
tion of the jet. Rather, the dominant coUimating force arises 
far from the axis between ^? ~ 5 — 15° where both B4, and Ip 
are well behaved. In fact, the gradients of (r sin OB^Y ^rid Bp 
both contribute to the magnetic force and in §4.3, we discuss 
in detail how these magnetic fields launch and coUimate the 
outflow. 



4 FLUXES, FORCES, AND VELOCITIES IN THE JET 
4.1 Matter and angular momentum flux 

4.1.1 Matter flux 

Here we analyze the matter and angular momentum fluxes in- 
tegrated over the star's surface (r = 1) and the r = 20 surface 
for the reference simulation. The mass accretion rate through 
a given surface is 

M = j dS-pvp = I dSpvr , (7) 

where dS is an outward facing surface area element and Vr 
is the radial component of the velocity. On the stellar sur- 

^ In contrast, the opposite is true at small radii (r < 5) and angular 
momentum is primarily transported away by the magnetic field. 



face, we measure the accretion rate Ms by only considering 
the inward (,Vr < 0) matter flux; conversely, for the outer r = 
20 surface we measure the outflow rate in the jet Mw by only 
considering the outward Vr > matter flux. Fig. 7a shows the 
inward Ms across the stellar surface and outward Mw across 
the r = 20 surface as a function of time. The total matter 
flux onto the star Ms slowly increases from the time the disk 
first reaches the star (t = 300) until it stabilizes around t = 
1300. Similarly, the total matter flux through the jet stabilizes 
slightly earlier at t = 1200. 

4.1.2 Angular Momentum Fluxes 

The angular momentum flux density consists of three com- 
ponents: angular momentum carried by the matter Fmi, by 
the magnetic field FLt, and by the viscosity Flv The total 
angular momentum flux density is therefore 

Fl = Flmi + FLt + Flv, 

where 

Fliii = r sin 6pv^Vp , 

Fl* ^-rsin^:^ , (8) 
47r 

Flv = -vtp{rs.mefVVt , 

with the last term non-zero only inside the disk. As with the 
matter flux, we integrate the angular momentum fluxes at the 
star (r = 1) and at the r = 20 surface, 

L = dS • Fl = JdS- (FLm + FLf + Flv). (9) 

Fig. 7b shows the integrated angular momentum flux as a 
function of time. and L^/ measure the angular momen- 
tum being added to the star by matter and magnetic fields. 
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Figure 7. Matter and angular momentum fluxes onto the star and into thie outflow. Panel (a) : Matter flux from the disk onto the star Ms and into 
the jet Mto as a function of time. The bottom plot shows the absolute value of the accretion-to-ejection ratio |Mt„/Ms| as a function of time. Panel 
(b): Angular momentum flux rates added onto the star by matter Lg,n and by the magnetic fields L^y . Also plotted is the angular momentum loss 
rate due to matter carried away by the jet and by the magnetic fields L^^f. The bottom plot shows the ratios of angular momentum fluxes 
leaving in the jet to the angular momentum fluxes onto the star. 



respectively. The accreting matter adds angular momentum 
directly (Lam) and spins up the star. The magnetic fields also 
allow the matter to spin up the star indirectly: the magneto- 
spheric radius Xm ~ 1-1 (calculated by equating the magnetic 
and gas pressure in the disk) is smaller than the corotation 
radius rcor = 3 and hence the inner disk drags the stellar 
dipole field, causing the star to spin up at a nearly constant 
rate (L^/). 

At the r = 20 surface, only the outward angular mo- 
mentum carried by the jet is calculated. The outward angular 
momentum transport due to magnetic fields, Lwf, is negligi- 
ble compared to the the angular momentum being carried 
away by the matter, Lwm. Like the matter flux, the angular 
momentum flux in matter becomes steady around t = 1200. 
Roughly a third of the disk's angular momentum exits the sys- 
tem through the jet, while only a small fraction of the incom- 
ing angular momentum is added to the star. The majority of 
the disk's angular momentum is transported by viscosity back 
and out through the disk itself. 



4.2 Velocities in the jet 

Fig. 5 shows a few of the velocity streamlines flowing from the 
base of the jet to the edge of the simulation region. The figure 
shows that the launching region of the jet is localized to the 
very inner disk. To probe the structure of the jet in detail, we 
analyze the velocities within the jet at t = 860. Fig. 8 shows 
velocities along a "representative" field line running through 
the jet in the t = 860 reference case. The field line originates 
from the disk at r = 1.24 and extends upward through the 
jet to the outer boundary of the simulation region (for a vi- 
sual representation, see the denoted field line in Fig. 9). We 
define the new coordinate s which traces the linear distance 




Distance along field line, s 

Figure 8. Velocities as a function of distance along a representative 
field line. The field line is anchored in the disk at r= 1.24 and extends 
through the jet to the outer boundary of the simulation region. A to- 
tal of eight velocities are plotted. Starting from the left: Usm - slow 
magnetosonic velocity, v^p - poloidal Alfven velocity, - velocity of 
matter tangential to the field line, v± - velocity of matter perpendic- 
ular to the field line, - fast magnetosonic velocity, - toroidal 
matter velocity along the field line, and lastly, Vcsc - the escape veloc- 
ity 

along the poloidal field line with s = at the base^. Fig. 8 
shows that very close to the star, — the matter velocity tan- 
gential to the field line — is small; however, ?;|| quickly rises 
above the slow magnetosonic Vbiu, poloidal Alfven VAp, and 
fast magnetosonic ufm speeds as the matter moves away from 
the star, indicating that the flow is matter dominated past s ^ 
13; most of the acceleration occurs before the matter crosses 



^ Since the representative field line is nearly vertical, s ^ r cos 9 
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the fast magnetosonic surface. Fitting a power law to the re- 
gion of strong accleration gives oc s°'^^; past s ~ 5, the 
acceleration mechanism weal<ens and the slope flattens such 
that V|| oc s"^*^. Near the disk, the toroidal component of the 
velocity, v^,, is close to the Keplerian velocity meaning that 
the matter in the jet initially corotates with the disk. How- 
ever, outside the disk, -y^ falls off as and exceeds 
at s w 12 due to the continued acceleration of the matter. 

Far away from the star, the tangential velocity is nearly 
double the escape velocity Vcsc and the matter in the jet easily 
escapes from the system. In contrast, the perpendicular veloc- 
ity v± in the jet is quite small, indicating that matter diffuses 
across the magnetic field very slowly. The small oscillations 
perpendicular to the field line may be an Alfven wave, but 
this remains to be investigated. 



4.3 Launching and collimation mechanisms 

4.3.1 Forces in the jet 

The force per unit mass tangent to a poloidal magnetic field 
line can be calculated by taking the dot product of the Eu- 
ler equation with the b unit vector which is parallel to the 

poloidal magnetic field line Bp. Following the derivation in 
Ustyugova et al. (1999), the total force tangent to Bp is 



fu 



= b • (fp fc + fc + fM 

1 dP 9$ 



(10) 



p OS OS rsm 



-sine+ b- [(V X B) X Bl. 

6 Attp 



Here, fp, fc, fc, fM are the pressure, gravitational, centrifugal 
and magnetic forces, respectively; O is the angle between the 
poloidal magnetic field line and the axisymmetry axis; $ is the 
gravitational potential; s is the previously defined coordinate 
which traces the distance along the poloidal field line. 

The pressure gradient force, /p y — —{1/ p)(dP/ds), 
dominates within the disk. The matter in the disk is approx- 
imately in Keplerian rotation such that the sum of the grav- 
itational and centrifugal forces roughly cancel (fc-HC ~ 0). 
Near the slowly rotating star, however, the matter is strongly 
coupled to the stellar magnetic field and the disk orbits at sub- 
Keplerian speeds, giving fc-i-c <, 0. The tangential magnetic 
force (the last term of Eq. 10) can be expanded as 



/i 



1 



b • [(V X B) X B] 



(11) 



47rp 

__ 1 d{rsmdB4,) 
87rp(rsin^)2 ds 

(Lovelace et al. 1991). Note that rsin^B^ is the normalized 
poloidal current flowing through a surface of radius r from 
colatitude zero to 9. 

The force per unit mass perpendicular to a poloidal field 
line is given by 



-I- 



-{Vp - "Ap)-^ 

iAfp+?i 

p dn \ Stt 

COS 



a(rsin6»B^)2 



87rp{r sin 6)^ 



dn 



9$ 

dn 



(12) 



r sin 6 

where n is a coordinate normal to the poloidal field line 
(Ustyugova et al. 1999). In §4.2 we showed that far from the 



disk, the poloidal velocity in the jet Vp dominates over the 
poloidal Alfen speed VAp (see Fig. 8), permitting us to ignore 
the poloidal Alfven term. Additionally, far outside the disk, 
the pressure gradient term is negligible and gravity is weak 
enough to be ignored. With these simplifications, Eq. 12 re- 
duces to 



f 2de 

hot,± - ~^p-Q^ 



1 dBl 
8np dn 
dirsineB^f 



V<4 COS 6 



8np{r sin 6)^ 



dn 



(13) 



Once the jet begins to collimate, the curvature term 
—VpdO/ds also becomes negligible. The magnetic force may 
act to either collimate or decoUimate the jet, depending on 
the relative magnitudes of the toroidal (rsin^i?^)^ gradi- 
ent (which colUmates the outflow) and poloidal Bp gradient 
(which "decoUimates"); in our simulations, the collimation 
of the matter implies that the magnetic hoop stress is larger 
than the poloidal field gradient. Thus the main perpendicular 
forces acting in the jet are the coUimating effect of the toroidal 
magnetic field and the decoUimating effect of the centrifugal 
force. 



4.3.2 Analysis of the forces 

Here, we analyze the role that the magnetic, gravitational, 
centrifugal, and pressure forces play in launching and coUi- 
mating the outflow. Fig. 9 shows vectors denoting the mag- 
nitude and direction of the various forces along a "represen- 
tative" field line running through the jet: this field line is the 
same as the one used previously for analysis of the velocities 
in §4.2 and is anchored in the disk at r = 1.24. Fig. 10 shows 
the tangential and perpendicular projections of the individual 
forces onto this field line. 

Launching mechanism The top row of Fig. 10 shows several 
different views of the tangential projections of the forces onto 
the representative field line. In order to accelerate the matter 
and launch an outflow, the sum of the tangential forces in the 
jet /tot, II must be positive overall. The tangential projections 
of the gravitational and centrifugal forces are nearly equal 
and opposite everywhere in the jet (i.e. /g-i-c,|| ~ 0) except 
near the star where gravity dominates and acts to hold the 
matter in the disk (e.g. Fig. lOIa). In contrast, the magnetic 
pressure force /m.|| acts in opposition to gravity and serves 
to accelerate the matter out of the disk along the magnetic 
field lines. Fig. lOIb shows that the gas pressure force /p_|| 
is small outside the disk and plays no role in launching the 
outflow. Since /g+c,|| and /p y are both negligible far from 
the star, the magnetic force dominates within the jet giving 

/tot, II ~ /m,|| • 

The effect of the various forces can be seen in the three 
panels of Fig. 9: fa+c pulls the matter downward toward the 
disk (Fig. 9b) while the the magnetic force /m acts to drive 
matter out from the disk (Fig. 9c). When these forces are 
summed together, the resulting force /tot shown in Figs. 9a 
and lOIb is positive overall; note that even far from the disk, 
/tot, II > and the matter is continually accelerated. Since 
the centrifugal force is completely canceled by gravity, the jet 
is driven by a purely magnetic force and hence the launch- 
ing mechanism is purely magnetic as well. This mechanism 
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Figure 9. Forces along a field line in the jet. Panel (a) shows the poloidal matter flux density p |vp| as a background overplotted with poloidal 
magnetic field lines. The vectors show the total force ftot along a representative field line originating from the disk at r = 1.24. Panel (b) plots 
the angular velocity U as the background. The vectors show the sum of the gravitational + centrifugal forces fo+c along the representative field 
line. Panel (c) shows the poloidal current Ip as the background. The vectors show the total magnetic force along the representative field line. 
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Figure 10. Force components along a field line in the jet. The top row of panels (la), (lb) & (Ic) show the tangential projection of the forces 
along the representative field line in: (la) the region just above the disk; (lb) the whole simulation region; (Ic) the region far above the disk. 
A positive tangential force indicates acceleration away from the equatorial axis (i.e. the disk plane). Panels (Ila), (lib), & (lie) show similar 
plots of the perpendicular projection of the forces along the representative field line. A positive normal force indicates acceleration away from the 
axisymmetry axis. 



is similar to the inner disk wind model discussed in Lovelace 
et al. (1991) and observed in simulations of conical winds in 
R09. 

Collimation mechanism As discussed previously in §4.3.1, 
the collimation of the jet is a competition between the decol- 



limating effect of the centrifugal force (see Fig. 9b) and the 
coUimating force of the magnetic hoop-stress (Fig. 9c) ; in or- 
der to coUimate the jet, the sum of these perpendicular forces 
must be negative overall. The bottom row of Fig. 10 shows 
the normal projection of the various forces onto the repre- 
sentative field line. Very close to the star, the perpendicular 
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Figure 1 1 . Dependence of the angle between the reference poloidal 
field line and the axisymmetry axis on the distance along the field 
line. 

component of the centrifugal force is large and the net force 
is positive, pushing some of the matter into a "spur" (see Fig. 
9a). However past s = 2.5, the centrifugal force weakens as 
the matter moves away from the rotation axis, allowing the 
magnetic hoop-stress to coUimate the matter into a jet. Fig. 
9 shows this very clearly: the centrifugal force in Fig. 9b acts 
opposite to the magnetic force in Fig. 9c. When the two forces 
are summed (Fig. 9a), the forces largely cancel and the resid- 
ual points inward toward the axis, serving to coUimate the 
outflow into a jet. Note that the scale of the vectors in Fig. 
9a is much smaller than the scale in in Figs. 9b and 9c: the 
collimation is a delicate balance between the magnetic and 
centrifugal forces in the jet. 

Fig. 11 shows the dependence of — the angle between 
the reference poloidal magnetic field line and the axisym- 
metry axis — on the distance along the field line s. For the 
poloidal field line in the reference case, we measure a half- 
opening angle 6 w 4° for s > 25. 



5 DISCUSSION 

5.1 Applications to young accreting stars 

The dimensionless form of the simulations permit us to ap- 
ply the results to a variety of stars of different sizes, from 
neutron stars to young T Tauri stars. In R09, applications of 
the conical winds to different tjrpes of stars were discussed 
in detail and more recently, a model similar to the one pre- 
sented in this paper has been applied to the rapidly accreting 
FU Orionis (Konigl et al. 2011). In our simulations, the high 
M disk compresses the magnetosphere almost to the surface 
of the star, and the model is well suited for understanding 
young stars undergoing periods of enhanced accretion such 
as FUORs, EXors and magnetized CTTSs. 

EXors EXors represent an interesting stage in the evolution of 
young stars where the accretion rate may be enhanced up to 
M ~ 10~® — 1O~^M0 /yr resulting in the ejection of powerful 



outflows (e.g. Coffey et al. 2004; Brittain et al. 2007). The ac- 
cretion rate obtained in our simulations is M = MMo (where 
M fa 160 is the dimensionless accretion rate obtained in sim- 
ulations, see Fig. 7). The reference accretion rate Mo (and the 
dimensional accretion rate) strongly depend on the magnetic 
field of the star. For EXors, we choose the value B» = 800 G 
and obtain M ~ 1.1 x lO~''M0/yr. Table Al shows the other 
reference values for this case. 

FUORs The model is also well suited for describing FUORs 
for which the accretion rate is very high, M 10^^ — 2 x 
lO^^MQ/yr, and the disk can strongly compress the magne- 
tosphere. In Konigl et al. (2011) the properties of the wind 
(density, velocities) were compared with those derived from 
the blue-shifted spectral features observed in the wind of FU 
Ori (Calvet et al. 1993; Hartmann & Calvet 1995) and reason- 
able agreement was found. Table Al shows an example of the 
reference values for FUORs. Here, we increased the magnetic 
field of the star up to S* = 2.5 x 10^ G in order for the magne- 
tosphere to be able to stop the disk near the star. For this B*, 

we obtain an accretion rate of M = MMo = 2.5 x IO'^Mq /yr 

6 

CTTS The model can be also appUed to young stars — such as 
protostars and CTTS — with much lower accretion rates. Com- 
pared to FUORs, CTTSs show relatively small accretion rates 
of M « 10"^ - lO"*M0/yr. Typical magnetic field strengths 
around T Tauri stars are on the order of one to a few kilo- 
gauss (see Yang & Johns-KruU 2011). However, the stellar 
field may be dominated by higher-order components of the 
magnetic field while the dipole field (which dominates in the 
disk) may be comparatively weak (Donati et al. 2007; Gre- 
gory et al. 2008; Donati et al. 2011). Taking this into consid- 
eration, we choose a weak equatorial dipole field of B,. = 100 
G as a CTTS reference value (see Table Al) and obtain a disk 
accretion rate of M = MMo = 1-8 x lO'^'Me/yr. 

The ejection-to-accretion ratio obtained in our model, 
Mwind ~ 0.2Maccr (sec Fig. 7), is somewhat larger than the 
typical values observed in CTTSs (<0.1) or FUORs (~0.1) 
(Hartmann & Kenyon 1996; Coffey et al. 2008). In our sim- 
ulations, the jets are successfully launched when the initial 
outflow is able to penetrate through the matter-dominated 
corona. The coronal density in our simulations is quite low 
(typically ~ 10^'' of pd). Despite this, the corona is initially 
matter-dominated because the dipole magnetic field strength 
decreases rapidly with distance. The many simulations which 
produced collimated jets (see Tab. 1) showed that the mag- 
netic launching mechanism is robust and can drive jets from 
disks with a wide range of densities and accretion rates. How- 
ever, only the most powerful outflows can penetrate through 
the matter-dominated corona. It is possible that we observe a 
selection effect in which only most powerful jets successfully 
penetrate through the dense corona, and in reality (when the 
coronal density is even lower), weaker jets can be launched. 
In that case, we expect the matter ejection-to-accretion ratio 
Mwind /Maccr to be lower. 



Note that Konigl et al. (2011) adopted a larger stellar radius and a 
smaller M, finding =2.1 kG. 
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5.2 Collimation in different simulations 

In Tab. 1 we show a sample of the simulations which pro- 
duced coUimated jets. The simulations exhibited a wide range 
of degrees of collimation with half-opening angles ranging be- 
tween 6 = 4° and 6 = 20° . The collimation observed in the 
jets is stronger than in the conical winds studied in R09^. 

The key difference between the simulations presented in 
R09 and the new simulations presented here is the size of the 
simulation region: the new simulation region is nearly triple 
the size of the regions in the previous work. The larger simu- 
lation region in our model setup results in a higher accretion 
rate*^ . In addition, in the larger region, incoming disk matter 
collects and compresses a larger portion of the star's magnetic 
flux toward the star. Both factors lead to a higher magnetic 
pressure at the disk-magnetosphere boundary and to a cor- 
respondingly higher magnetic force. This results in the emer- 
gence of a more powerful, more magnetized outflow, lead- 
ing to stronger collimation of the jet. We measure magnetiza- 
tion level with the kinetic plasma parameter I3i which takes 
the ratio of the thermal plus ram pressure to the magnetic 
pressure (Eq. 6). In all of the cases shown in Tab. 1, the jets 
are matter-dominated (i.e. /3i,jct > 1). However, in strongly 
matter-dominated jets (/3i,jet > 10), weaker collimation is 
observed — for example, in the conical winds observed in R09, 
O ~30° and /3ijet = 10 — 30. In contrast, when the outflow 
is only slightly matter-dominated (like in the reference case 
of this paper where /Jijet « 1 — 6, 8 « 4°) then collimation 
is strong. In §3.2 we argue that the spurious toroidal field at 
6 < 3.6° has no influence on the jet collimation. 

Recently, Kurosawa et al. (2011) performed a special set 
of simulations for a magnetized star with much smaller ac- 
cretion rates (as in a CTTS) where the disk does not com- 
press the magnetosphere as strongly and is truncated at a few 
stellar radii. These simulations are much more applicable to 
CTTSs which usually have large magnetospheres and showed 
that outflows are usually less coUimated when the accretion 
rate is low. 



5.3 Comparison with observation 

Recent observations of several CTTSs with outflows have 
shown that ~10 AU away from the star, the outflows have al- 
ready become coUimated. However, in almost all of the cases 
the degree of collimation is low, with opening angles of ~20- 
30° at distances 10-50 AU (e.g. Hartigan et al. 2004). These 
outflows eventually become very well coUimated at larger dis- 
tances, with opening angles of 2 — 5° further out (Dougados 
et al. 2000). The jets in our simulations exhibit much faster 
collimation — 4 — 20° at 0.4 AU and are likely most relevant 
to highly accreting magnetized stars such as FUORs or EX- 
ors where the outflow is expected to be launched locally from 
a highly compressed disk-magnetosphere boundary. However, 
observational coverage of jets around FUORs and EXors is still 



^ We should note that even in the case of the conical winds, the half- 
opening angle at the boundary reaches 20° — 25° at the end of the 
simulation run (see Fig. 2, t = 700 panel) which is marginally close 
to the least coUimated jets observed in current simulations. 
* At the external boundary we fix the density of the incoming matter 
and hence M increases with the size of the simulation region. 



sparse and there exist few constraints on the sub-AU struc- 
ture of outflows from these highly accreting stars. However, as 
mentioned previously in §5.1, comparisons of the our model 
with spectrally-derived properties of the wind in FU Ori star 
show that the simulations do agree with observed outflow ve- 
locities (Konigl et al. 2011). 

Alencar et al. (2005) performed a spectral analysis of the 
Ha, Hp, and NaD emission in the magnetosphere and at the 
base of the outflow in the highly accreting CTTS RW Aur A (M 
~ lO"'^-® to lO"*^ Mo/yr, Hartigan et al. 1995; White & Ghez 
2001). Their results suggest that the high- velocity microjet in 
RW Aur A is likely launched as a "narrow wind" from a disk 
region between 2.2-5 R, around the star, consistent with the 
launching region observed in our simulations. Their spectral 
fits also show that the field lines at the base of the outflow are 
inclined at 6 = 30 — 40°, consistent with the jet launching 
angle in our models (see Fig 11). 

High velocity outflows launched from the inner disk 
(such as the coUimated jets in our simulations) may explain 
the central component of jets with a layered "onion-skin" 
structure where the highest velocity, well coUimated portion 
of the outflow flows close to the axis (Bacciotti et al. 2000) . 
If the disk is also threaded by an ordered magnetic field, then 
an extended disk wind may be responsible for the outer, less- 
coUimated, low- velocity layers (Ouyed & Pudritz 1997; Fer- 
reira et al. 2006; Fendt & Elstner 2000). However, the sim- 
ulations show a broad range of collimation and it is possible 
that the more weakly coUimated jets are related to these outer 
layers as well. 



6 CONCLUSIONS 

Through 2.5D MHD simulations, we have achieved robust, 
fully coUimated jets emerging from the disk-magnetosphere 
boundary of an accreting magnetized star. We observe purely 
magnetic launching and collimation mechanisms: the cold 
matter at the disk-magnetosphere boundary is driven up- 
ward by the magnetic pressure and coUimated by the heli- 
cally wound magnetic field lines extending up from the disk. 
Approximately one-fifth of the incoming disk matter exits the 
system through the jet; the jet also advects some angular mo- 
mentum out of the disk but the majority is transported out- 
ward through the disk by viscous stress. The degree of colli- 
mation in the jet is likely connected to the level of magneti- 
zation within the jet: our simulations show that strongly mat- 
ter dominated jets with fli > 10 are less coUimated than the 
weakly matter dominated outflows where /3i > 1. In the refer- 
ence simulation presented in this work, the matter in the jet is 
only slightly matter-dominated (/3i J> 1) resulting in a mag- 
netic field that strongly coUimates the jet to a half-opening 
angle of « 4° at the top of the simulation region. Previous 
simulations by our group performed in a smaller simulation 
region showed the emergence of weakly coUimated conical 
winds with a larger jii in the outflow (R09) . The high levels 
of magnetization necessary to produce a strongly coUimated 
jets may arise around stars such as EXors or FUORs where a 
high accretion rate strongly compresses the magnetosphere. 

Observations suggest that the jets from protostellar sys- 
tems become coUimated at distances less than 10 AU away 
from the star (the present resolution limit) (Hartigan et al. 
2004; Coffey et al. 2008). Our simulations show that the jet 
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is launched and collimated on scales which are on the order 
of tens of stellar radii. However, the jets presented here may 
represent just the well collimated, high-velocity core of the 
"onion-skin" type outflows observed around young stars. 
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CTTS 


EXors 


FUOri 


initial 


M*[Mq] 


0.8 


0.8 


0.5 


R* [R©] 


2 


2 


2.5 


U [days] 


1.90 


1.90 


3.37 


B* [G] 


100 


800 


2500 


derived 


Ro [cm] 


1.39 X 10" 


1.39 X 10" 


1.74 X 10" 


to [days] 


0.366 


0.366 


0.648 


vo [cms~i] 


2.76 X 10'^ 


2.76 X lO'' 


1.95 X lO'^ 


Bo [G] 


10 


100 


250 


fo [dynes/gm] 


5.48 X 10^ 


5.48 X 10^ 


2.19 X 103 


7o [G cmVs] 


2.09 X 10^2 


1.67 X 10^3 


6.52 X 10^3 


/() [A] 


6.96 X 10^2 


5.57 X 10" 


2.18 X 10" 


Po [g cm-'^] 


1.31 X 10-13 


8.39 X 10-" 


1.64 X 10-1'' 


no [cm^'*] 


7.83 X 10^0 


5.01 X 10^2 


9.79 X 10^3 


Mo [Mq yr-i] 


1.11 X 10"'' 


7.12 X 10-* 


1.54 X 10-** 


Eo [ergs-i] 


5.35 X 10^1 


3.43 X 10-" 


3.70 X 103* 


Lo [ergs-i] 


2.70 X 1035 


1.73 X 1037 


3.29 X 1038 


Td [K] 


9178 


9178 


4589 


Tc [K] 


9.18 X 10^ 


9.18 X 10® 


4.59 X 10« 



Table Al. Calculated reference values for various types of young, ac- 
creting stars. We choose the values of the stellar mass M«, radius R*, 
period t*, and equatorial magnetic field B, of the star and derive 
the other reference values from these initial values. The stellar rota- 
tion period is set by the corotation radius of the disk which is set to 
Rcoi = 3. The magnetic field B, is chosen such that the disk stops 
close to the stellar surface. To apply the simulation results to a partic- 
ular class of stars, one needs to multiply the dimensionless values by 
the reference values shown in this table. 



Yang H., Jolins-Krull C. M., 2011, ApJ, 729, 83 
Zanni C, Ferrari A., Rosner R., Bodo G., Massaglia S., 2007, 
A&A, 469, 811 



APPENDIX A: REFERENCE UNITS 

We take the reference mass Mo to be tlie mass M, of tlie 
star. The reference radius is taken to be the radius of the star, 
Ro — R* ■ The magnetic field is a dipole field and its value at 
the equator, B,, is chosen such that the accretion disk (with 
typical accretion rates for young stars) can compress the mag- 
netosphere almost to the stellar surface. The reference mag- 
netic field is Bo ~ B* //i, where /i is the dimensionless mag- 
netic moment. The reference velocity is vo = {GM / RqY^'^ , 
the Keplerian velocity at the surface of the star. We measure 
time in units of to = 27rJ?o /vo (the Keplerian rotation period 
at r = Ro). The reference force per unit mass is vo/Ro- The 
reference density is taken to be po = Bq/vo- The reference 
temperature is To = Po/ijtpo) = vo/TZ, where TZ is the gas 
constant. The reference mass accretion rate is Mo = povoR^. 
The reference energy flux is Eo = MoVq- The reference angu- 
lar momentum flux is Lo — MovoRo- The reference current is 
lo = cRoBq. Tab. Al shows examples of reference values for 
several different types of young accreting stars. 



APPENDIX B: VISCOSITY AND DIFFUSIVITY IN THE DISK 

The momentum flux-density tensor T in the MHD equations 
(Eq. 2) can be expanded as 



Tik = pviVk + PSik + ( ~ 



+ Tik. 



(Bl) 



Here, P is the gas pressure, 5ik is the Kronecker delta func- 
tion, and Tife describes the viscous stress from small-scale tur- 
bulence in the velocity and magnetic fields. This non-ideal 
portion of the stress tensor, nk, is dominated by two compo- 
nents 

Tr4, = —vtprsmO — , res = —ftpsmO-p— , (B2) 
or or 

where Q, is the angular velocity and vt is the coefficient of 
the kinematic turbulent viscosity. In our simulations, we do 
not consider the effects of viscous heating or radiative cooling 
based on the expectation that the two processes will compen- 
sate for one another. Rather, the main role of viscosity is to fa- 
cilitate the transport of angular momentum outward through 
the disk, thereby permitting matter to accrete inward toward 
the star. 

To estimate the turbulent viscosity vt, we adopt the 

Shakura & Sunyaev (1973) a model which approximates the 
viscosity coefficient as vt = av(?s£,/^K, where Cs — {P j pY^"^ 
is the isothermal sound speed, ^ is a density threshold coef- 
ficient, and f2/f is the Keplerian angular velocity. Similarly, 
we estimate the turbulent diffusivity coefficient in Eq. 3 as 
r]t = oid(?g^l^K (Bisnovatyi-Kogan & Ruzmaikin 1976). Both 
a„ and ad are dimensionless coefficients which are treated as 
parameters in our model. ^ is a coefficient which sets a density 
threshold for the viscosity and diffusivity 







PA 
1 



if p < pd/4 
if pd/4 < p < Pd 
if p > Pd 



-Pd/4 



(B3) 



where pa is the disk density parameter at the external bound- 
ary (see Tab. 1). The density threshold coefficient ^ varies 
between to 1 and acts to smoothly "turn on" viscosity and 

diffusivity for regions with p > pd/4. Fig. Bl shows the vis- 
cosity and diffusivity profiles of the inner disk of the reference 
simulation at t= 860. Recall that the reference case used coef- 
ficents of ad = 0.1 and a„ = 0.3 for the turbulent diffusivity 
and viscosity. At each timestep, the turbulent diffusivity and 
viscosity profiles are evolved and reapplied to every point on 
the grid to ensure self-consistency of the simulation. Above 
the Pd /4 (^ > 0) threshold, viscosity becomes important and 
we include rit and vt by numerically integrating the non-ideal 
MHD equations using an explicit conservative Godunov-type 
numerical scheme. Outside the disk where the density is low, 
5 = and we solve the equations of ideal MHD. 



APPENDIX C: CURRENT FLOW IN THE SIMULATION 
REGION 

It is important understand the nature of the current flows in 
the simulation region. For general axisymmetric MHD flows 
the poloidal magnetic field (,[Br, Be] in spherical coordinates) 
can be written as 



Bp = V X 



e 

rsinO dr ' 



(CI) 
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Figure Bl. Zoomed in plot of the diffusivity and viscosity profiles near 
the star at t = 860. In the 5 = region, the turbulent viscosity ut and 
diffusivity rjt are turned-off. Above the Pd /4 density threshold, both 
are smoothly turned on as § — > 1 and we solve the non-ideal MHD 
equations in the disk. 

where *(r, — r sin 9 Atf,{r, 6) is the "flux function" for the 
poloidal magnetic field and is the vector potential (e.g., 
Lovelace et al. 1986). The "i/ = const lines label the poloidal 
field in that Bp ■ V* = 0. 

For axisymmetric and non-relativistic MHD flowrs, 
Ampere's law gives the poloidal current density (in cgs units) 
as 

J. = f V X iB,4>) -^'4^—/^^ (C2) 
47r r'' sin oO r sm or 

where H{r,6) = cr sin 6 Brf,{r, 9)/ Atv acts as a the "flux func- 
tion" for the poloidal current density in that Jp • \7H = (e.g., 
Lovelace et al. 1986). 

In cgs units, the total poloidal current through a spherical 
cap, r = const and colatitude from to 9, is simply 

/■» 

^p.cgs = 2tv r sin9d9Jpr = 2iTH{r,9) — -r sin9B^(r,9), 
Jo 2 

(C3) 

where Jpr = f • Jp. In all the plots we show the normalized 
current, /p = r sin9B^{r,9), which is proportional to /p,cgs. 
For the dipole-type field symmetry about the equatorial plane 
assumed in our simulations, ^'(r, 9) = '^{r,-K — 9); that is, '3/ 
is an even function about the equatorial plane. This requires 
that the toroidal magnetic field be an odd function, B^{r, 9) — 
-B4,{r,-K- 9) so that B^{r,n/2) = (Lovelace et al. 1987). 
Thus we necessarily have 

7p(r,7r/2) = 0. (C4) 

Thus, our MHD simulations guarantee that the net current 
flow through the upper (or lower) hemisphere is exactly zero. 
A current outflow in a jet with, for example, ^ < 1 is exactly 
balanced by a current inflow in the region tv/2 — 9 < 1. How- 
ever, in the general case of no symmetry about the equatorial 
plane, we may have Ip{r, tv/2) / (see Lovelace et al. 2010). 

The top panel of Fig. CI shows the poloidal current, 
Ip{r,9), surfaces inside the simulation region. The values of 
Ip in normalized units are shown on the lines. Fig. Clb shows 




Figure CI. Panel (a): The color background shows the values of 
Ip{r, 9) and the lines are current density lines. The numbers on the 
contours denote the local value of the current in normalized units. 
Panel (b) : The 9 cross-section of the current and /3i at r = 20, (see 
dashed line in the top panel) . The dotted vertical lines represent the 
/3i = 1 surfaces. 

that the largest current flow is between the jet and the disk 
with a value Ip — 3, corresponding to 2.1 x 10^^ Amperes for 
a disk around a CTTS (see Tab. Al). The current flow inside 
the jet is almost parallel to the jet, and is smaller than the 
current within the disk. The current carried by the jet flows 
out through the outer spherical boundary. At the same time, 
an equal amount of "return" current flows inward through 
the outer boundary and along the surface of the disk. The jet 
current and its return current have opposite signs and they 
repel as a result of their magnetic interaction mediated by 
the toroidal magnetic field B,j,. This repulsion between the jet 
and its return current has been clearly observed and discussed 
in previous MHD simulations by Ustyugova et al. (2000) and 
Nakamura et al. (2008). The toroidal magnetic field is respon- 
sible for collimating the jet as discussed in §4.3. At the same 
time the toroidal field gives an outward radial force on the 
annular return current pushing it away from the jet axis. 
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